Observation of subdiffusion of a disordered interacting system 
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We study the transport dynamics of matter-waves in the presence of disorder and nonlinearity. 
An atomic Bose-Einstein condensate that is localized in a quasiperiodic lattice in the absence of 
atom-atom interaction shows instead a slow expansion with a subdiffusive behavior when a con- 
trolled repulsive interaction is added. The measured features of the subdiffusion are compared to 
numerical simulations and a heuristic model. The observations confirm the nature of subdiffusion as 
interaction-assisted hopping between localized states and highlight a role of the spatial correlation 
of the disorder. 
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The combination of disorder and nonlinearities deter- 
mines the transport properties of many physical systems, 
including metals and superconductors [l|, graphene (3 
and DNA or light in disordered nonlinear media [J- 
8] and ultracold quantum gases While a full under- 
standing of the interplay of disorder and nonlinearities 
has long been sought, systematic experimental investi- 
gations arc difficult and various aspects of this interplay 
have not been fully clarified. One of the open questions 
concerns the dynamics of a wavepacket expanding in a 
disordered potential in presence of a nonlinearity. There 
have been various predictions in theory and debate based 
on the results of numerical experiments over the last 20 
years about this subject [343, 0, El, E, [T3-IH, El] ■ Most 
authors agree that the nonlinearity should prevent local- 
ization and that the wavepacket should expand in a way 
which is slower than normal diffusion. However, experi- 
mental evidence of such subdiffusive expansion is to date 
still missing. 

Here we study the dynamics of Bose-Einstein conden- 
sates with controllable nonlinearity expanding along a 
one-dimensional quasiperiodic lattice. Despite its large 
spatial correlation, this kind of potential is known 2(], 21 1 
to feature exponentially localized states that are equiv- 
alent to those appearing in lattices with uncorrelated 
disorder described by the Anderson model [22j. Such 
a system has been successfully exploited for the investi- 
gation of the delocalizing effect of repulsive interactions 



on a trapped system in equilibrium [23], |24j. If a non- 
interacting gas is let free to expand along the quasiperi- 
odic lattice, no transport is observed |25j because all 
single-particle eigenstates are localized. By adding a con- 
trolled interatomic repulsion, we now observe a slow in- 
crease of the width a of the sample that asymptotically 
follows a subdiffusive law: a(t) cx t a , with a — 0.2 — 0.4. 
We find that the exponent increases with the interac- 



tion energy, in qualitative agreement with both numerical 
simulations based on a ID discrete nonlinear Schrodinger 
equation (DNLSE) of the quasiperiodic lattice Q and the 
predictions of a heuristic model. Our observation con- 
firms the nature of subdiffusion as interaction-assisted 
hopping between localized states. The observed expo- 
nents are however larger than the one calculated for un- 
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correlated disordered potentials [1H3, 
suggesting a role of the spatial correlation of the disorder. 

The one-dimensional quasiperiodic potential is created 
by perturbing a primary optical lattice with a weaker 
incommensurate lattice [26j: V{x)—V\ cos 2 (kix) + 
V2 cos 2 (fc2x). Here ki—2-K/Xi are the wavevectors of the 
lattices (Ai=1064.4 nm and A2=859.6 nm). This poten- 
tial is characterized by the spacing d—Xi/2 and the tun- 
neling energy J of the primary lattice, and by the dis- 
order strength A, which scales linearly with V2 [27j . In 
the case of non-interacting particles this system consti- 
tutes an experimental realization of the Harper or Aubry- 



Andre model [20|, |2l| which shows a transition between 



extended and localized states for a finite value of the 
disorder A/ J = 2 [il l27j|. Above this threshold all 



single-particle eigenstates of the first band of the lat- 
tice are exponentially localized, with a localization length 
fwd/ ln(A/2J). 

We employ a Bose-Einstein condensate of 39 K atoms 
in their lowest internal state, whose s-wave scattering 
leng th a can be tuned by means of a Feshbach resonance 
2Sl 29|. The condensate is produced in an optical trap 



at a = 280ao, and contains about 5 x 10 4 atoms. We 
first load the interacting condensate into a quasiperiodic 
lattice with a constant A s» 3 J. The trap potential has 
a radial (axial) frequency of 27r x 50 (70) Hz, while the 
lattice beams give an additional radial confinement with 
frequency u r = 2n x 50 Hz. At a given time t = the 
optical trap is suddenly switched off, allowing the sam- 
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pie to expand along the lattice; at the same time, A and 
a are tuned to their final values within 10 ms, and kept 
there for the rest of the evolution. The subsequent change 
of the radially- integrated spatial distribution n(x) of the 
sample is then monitored by in-situ absorption imaging 
for increasing times, up to t = 10 s. The width of the dis- 
tribution is measured as the square root of its second mo- 
ment: <J={§ x 2 n(x)dx) 1 / 2 . The initial interaction energy 
per particle is estimated as E int —gN/2 J ip 4 d 3 x, where 
g=4:Trh 2 a/m is the coupling constant, tp is a Gaussian 
approximation of the single-site wavefunction and N is 
the mean atom number per site. We estimate that the 
initial distribution occupies on average 20±7 sites; this 
uncertainty translates into a 35% uncertainty on Ei nt . 
Note how in this work we stay in the regime where £ ~ d 
is much smaller than a. 
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FIG. 1: (color online) Time evolution of the width a for 
different initial interaction energies: Eint = (squares), 
Ei n t = 1.8J (triangles), and E in t = 2.3J (circles). The 
continuous lines are the fit with Eq.(l). The dashed lines 
show the fitted asymptotic behavior, while the dash-dotted 
line shows the expected behavior for normal diffusion. The 
lattice parameters are J/h — 180 Hz, A/J = 4.9. 

In Fig. [3] we present the typical time evolution of the 
width of the system for A/J > 2 for different interac- 
tion energies. In the absence of interaction, one observes 
only an extremely slow expansion, presumably due to 
technical noise. The introduction of a repulsive interac- 
tion allows the wave-packet to expand significantly: the 
expansion is however not ballistic since its velocity de- 
creases with increasing widths. To model the expansion 
we fit the measured evolution with the solution of a gen- 
eralized diffusion equation: 



a(t) = a (l + t/t ) a 



(1) 



which is expected to correctly model the overall behav- 
ior from short times to the asymptotic regime [3(J. Here 
do is the initial width, to is an "activation time" and a 
is the diffusion exponent (see also Eq.(2) later on). For 
normal diffusion a=0.5; from the fits instead we extract 
a w 0.2 — 0.4, indicating subdiffusion. From the fit we 



typically obtain values for to of the order of 1 s, indi- 
cating that the asymptotic behavior t > to is achieved 
approximately for one decade. The expansion becomes 
faster when increasing the interaction energy. This is 
confirmed by the systematic study of the exponent re- 
ported in Fig. [51 where one sees a clear increase of a with 
Eint, up to about a — 0.4. We compare these observa- 
tions with numerical simulations we have performed on 
a ID DNLS E Ml , for parameters close to the ones of the 
experiment [30j . The results shown in Figj6] are in qual- 
itative agreement with the experiment: a increases with 
Ei n t, to approach a saturation value around a ss 0.35. 
Note that noise in the experiment gives an a > already 
for Ei nt =0. This bias might justify the slightly larger 
values of a in the experiment than in theory j3(| • There 
is however a disagreement on the typical energy at which 
the saturation regime for the exponents is reached in the 
experiment and in the theory, which we will discuss later. 




FIG. 2: (color online). Diffusion exponent a vs the ini- 
tial interaction energy Eint in the experiment (triangles and 
squares) and simulations (circles). The experimental data are 
for A/J = 5.3(4) and two different values of the tunneling: 
J/h = 180 Hz (triangles) and J/h = 300 Hz (squares). The 
vertical bars are the fitting error of Eq.(l) to the data, while 
the horizontal bars indicate the statistical error. 

The observed subdiffusion confirms the microscopic 
mechanism of the expansion expected for an interact- 
ing disordered system where all the single-particle states 
are localized. The interaction breaks the orthogonality 
of the states and allows the transfer of population be- 
tween neighboring states. Since the transfer rate depends 
on Ei ri t , the velocity of expansion decreases as the sam- 
ple expands and becomes less dense. Various authors 
have built heuristic models for the expansion that re- 
late the transfer rate T to the variation of the width a 
in a variety of ID random models, such as the kicked 
rotor IlL the Anderson model and Klein-Gordon chains 

us mm. it is indeed possible to derive from basic 

principles, e.g. from perturbation theory, that the rate is 
proportional to the square of the interaction matrix ele- 
ments, r cx |F'| 2 , with \V'\ » gN/2(\I m \) dEij]. Here 
Iijki is an overlap integral between four single-particle 
states. The average (...) is extended to all quartets of 



3 



states that satisfy energy conservation, which for a dis- 
ordered system amounts to have the interaction coupling 
term \V[- kl \ — gN/2\Iijki\ larger than the energy separa- 
tion between the four states @, Hi- Note that \V^ kl \ is 
essentially Ei nt times a spatial overlap coefficient that is 
of order unity for neighboring states. 

The scaling of the atom number per site N with the 
width, N oc er -1 , leads to an equation of motion for a. 
For example, the idea followed for uncorrelated disorder 
[H, IH is to assume that T is essentially the instantaneous 
diffusion constant, i.e. da 2 /dt oc T. The resulting equa- 
tion da/dt oc a~ 3 can be integrated with initial condi- 
tions cr(0) = uq to get 



o-(i) =o-o(l + Ct/{aal /c 



(2) 



with a = 1/4 and C oc g 2 . This result holds in the regime 
of large interaction, where all the single-particle states 
within a few £ can be coupled by the interaction. Oth- 
erwise, also the number of possible couplings decreases 
with some power of gN and the expected exponent is 
smaller, e.g. a — 1/6 for uncorrelated disorder [3|. 

The values of a we extract from the experimental data 
and from numerical simulations are as large as 0.35-0.4, 
hence definitely larger than the maximum a = 0.25 found 
for the uncorrelated random models above. This suggests 
a role of the non-decaying correlation of the quasiperiodic 
lattice, originating from the quasiperiodic arrangement of 
the single-particle energies [27| . This might allow the ex- 
pansion to proceed through almost coherent hopping pro- 
cesses, in contrast to the incoherent mechanism expected 
for the uncorrelated models. In this case one might relate 
r to the instantaneous velocity, i.e. da/dt oc T, which 
leads to Eq. (2) with an exponent a — 1/3 in case of large 
interaction. This regime is reached when Ei nt is larger 
than the standard deviation of the single-particle ener- 
gies 5E. For Ei n t < SE, a numerical evaluation of the 
number of active couplings as in Ref . [2| shows that there 
is an additional scaling of V approximately as gN [3r3 |. 
Therefore in this regime one expects a smaller exponent 
around a = 1/4. This is consistent with what we obtain 
both in the experiment and in the simulations. 

Since the coupling between localized states is largest 
at the center of the sample, where the atom number 
per site is largest, one expects a faster expansion at the 
center than in the tails [1[. This behaviour is actually 
clearly visible in the typical shape of the cloud at a long 
t shown in Fig. [7K, which features an extended flat top 
and rapidly decaying tails. The flat-top shape gradu- 
ally emerges from initially gaussian profiles during the 
expansion, as can be seen in the typical evolution of the 
kurtosis, 7 = (J x A n(x)dx) / a 4 — 3, in Fig. [TJd. A com- 
parison with general models of nonlinear diffusion 0, [l6| 
is due in future studies. 

An additional feature of the experimental system is the 
presence of the radial degrees of freedom, something that 
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FIG. 3: (color online) a) Shape of an interacting sample after 
10 s of evolution; the dashed curve is a gaussian fit of the tails, 
b) Time evolution of the kurtosis 7; the dashed line marks the 
kurtosis of a gaussian distribution. The interaction energy is 
Eint ~ 1.8J, all other parameters as in FigEJ 



is absent in all numerical simulations and models devel- 
oped so far. In our case the radial quantum Tiui r R3 J/5 
is about one order of magnitude smaller than the typi- 
cal 8E w 3J. Therefore several radial states might be 
populated in presence of the interaction, which couples 
the various degrees of freedom. The radial excitation can 
be taken into account by inclu ding the radial wavefunc- 
tions in the overlap integrals / [30j ■ An increasing radial 
excitation can change the microscopic hopping dynamics 
in two ways: i) It increases the radial width of the sam- 
ple, and therefore decreases /, with a consequent slowing 
down of the expansion and a possible reduction of a. 
ii) It provides a smaller energy scale for the coupling of 
the localized states, since the axial energy mismatch can 
now be reduced in presence of a radial excitation. In- 
deed, while the minimum energy mismatch in ID is of 
the order of SE, it becomes of the order of hui r for large 
temperature samples in 3D. This effect tends to increase 
the overlap /, especially in the case of a small Ei nt that 
would not be able to couple efficiently states in ID. 

To check the relative magnitude of these two compet- 
ing effects we have compared the evolution of interacting 
samples at low radial temperature T r with those of sam- 
ples with the same g, but with a larger T r prepared via 
controlled parametric heating. The typical observation is 
reported in Fig. [U the expansion of the high-T r sample 
is faster, despite a reduction of Ei nt by approximately 
a factor 2, indicating that the second effect dominates. 
Raising the temperature therefore helps the coupling of 
a larger number of states by the interaction. In order to 
separate out the effect of the reduced Ei nt at the larger 
T r , we have also measured the time evolution of the sam- 
ple at low T r in which a similar Ei nt was achieved by re- 
ducing g. This measurement shows a slower expansion, as 
expected since now the reduced Ei nt is not compensated 
by an increased coupling efficiency. The effect of the ra- 
dial degrees of freedom can also explain the disagreement 
in the energy scales between theory and experiment in 
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Fig|6] in the ID simulation the relevant energy for being 
in the large interaction regime is SE w 3 J, while in the 
experiment in 3D this is shifted towards Tiui r m J/5. 

From a study of the radial momentum distribution of 
the interacting samples we also detect a radial heating 
during the 10 s evolution time, with a typical rate of 
the order of 3-5 nK/s that is at least twice as large as 
the heating of non-interacting samples [30j | . This heating 
leads to final T r that can be larger than SE/ks ~ 50 nK, 
so that it cannot be entirely justified by thermalization 
of the axial and radial degrees of freedom. We attribute 
such excess heating to energy transfer from the axial noise 
mentioned above to the radial degrees of freedom, medi- 
ated by the interaction. Finally, we note that for non- 
interacting samples we do not observe any dependence of 
the axial dynamics on the radial temperature, confirming 
the secondary role of the radial modes in the subdiffusion. 
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FIG. 4: (color online) Effect of a finite radial temperature. 
The scattering length and mean temperature of the three 
data-sets are: a = lOOao, T r = 60nK (circles); a — 400ao, 
T r = 60nK (triangles); a = 400a , T r = 200nK (squares). 
The continuous lines are the fit with Eq.(l), while the dashed 
lines are the fit to the asymptotic behavior. The fitted 
exponents are: a=0.32(2) (circles); a=0.34(2) (triangles); 
a=0.40(5) (squares). The lattice parameters are J//i=290 
Hz, A/J=3.9. 



In conclusion, our study gives evidence of the subdiffu- 
sive character of the expansion of an interacting system 
in a disordered potential. On the timescales considered 
here, it qualitatively validates the existing models of the 
subdiffusion as interaction-assisted hopping between lo- 
calized states, although it leaves open the question about 
the fate of the expansion at very long times [3, [35j . Our 
study raises new questions about the role of the corre- 
lations of the potential, which would be interesting to 
explore also in different kinds of disorder with intermedi- 
ate correlations. Also, we could drive our system into the 
quasi-lD regime by increasing the radial confinement to 
freeze out the radial degrees of freedom. There one could 
study the predicted temperature-induced metal-insulator 
transition [4] , or strongly correlated regimes 33|, |34 j . 
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SUPPLEMENTARY MATERIAL 
A PERTURB ATIVE MODEL OF SUBDIFFUSION 

In this section we discuss with some more detail the 
process of hopping between localized states induced by 
the interaction that gives rise to subdiffusion. A model- 
ing of the subdiffusion of disordered, interacting systems 
has already been given by several authors |l|-[3| by analyz- 
ing the properties of the Discrete Nonlinear Schrodinger 
Equation (DNLSE). Here we present an analogous heuris- 
tic model based on perturbation theory concepts. Al- 
though this model cannot be used to describe quantita- 
tively the full dynamics, it provides correct estimates for 
the subdiffusion exponents. Let us start by considering 
the localized eigenstates of the ID single-particle prob- 
lem, (pi, labeled by the site index i. These states have 
individual energies q with typical separation SE A 
and a characteristic localization length £ ps d/ln(A/2J). 
The interaction Hamiltonian operator is 




Axial position (arb. units) 

FIG. 5: Cartoon of the energy spectrum of the axial states 
along one quasiperiod of the quasiperiodic lattice (continu- 
ous lines) together with their associated radial states (dashed 
lines). The localization length £ is of the order of the lattice 
spacing d. 



W = I / 4>\x)ip f (x)4>(x)il>(x)dx , (3) 



where 



$( x ) = ^2b i ip l (x) 



(4) 



with bi the bosonic destruction operator at site i. In gen- 
eral, up to four states can be coupled by the interaction. 
The dominant process is the one in which two particles 
are moved from the states i and j with occupation N to 
two initially empty states k and I. Let us indicate by N 
the mean number of atoms per site in the initial distri- 
bution, which is related to the width in units of lattice 



sites a/d by the relation N — N to td/<7- The off-diagonal 
term of the interaction Hamiltonian for this process has 
the form 
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N / ip*(x)tp%(x)tpj(x)<pi(x)dx . 



(5) 



From perturbation theory we know that this process can 
happen at asymptotic times only if the energy is con- 
served, i.e. AE = | e/ + ej — €/■ — ej| = 0. In a disordered 
system it is impossible to have perfect energy matching 
of the states, but on the other hand if the transition rate 
for the process is large enough, even a finite energy differ- 
ence cannot be resolved and the transition can take place. 
As discussed for example in Refs. [3, [B[ the appropriate 
energy conservation requirement is actually 



\Vl ]U \>AE. 



(6) 



The transfer rate associated to this microscopic process 
is 



ijkl 



2?r \V- Jkl \ 2 
h AE 



(7) 



From Eq. 3 one sees that the coupling term is 
essentially the interaction energy Ei nt times an overlap 
integral I « exp(— L/£), where L is the mean separa- 
tion between the four states. We can therefore conclude 
that the macroscopic expansion of the system will be de- 
termined by all the microscopic processes between states 
that are laying within a few localization lengths from 
each other. It is instructive to numerically calculate the 
typical coupling strength in the quasiperiodic lattice for 
transitions between doublets of states, \V( ikk \/ AEik, for 
various separations between the two states. As shown 
in Fig. [5J the largest coupling is achieved on average 
for states that are separated by one quasiperiod of the 
lattice, d/(k2/ki — I). This result is a consequence of 
the peculiar spatial correlation of the eigenstates in the 
quasiperiodic lattice. 

One can distinguish two different regimes of expansion, 
along the lines of the discussion in Ref. Q and references 
therein. 

1) If the initial interaction energy is sufficiently large, 
then at least the strongest couplings of each localized 
state to the other states within one localization length 
are active. This regime is reached when Ei nt > SE. Here 
one can define a macroscopic rate T = (T^ki ) , where the 
symbol (...) denotes an average over the system exten- 
sion. From Eqs. 3 and 5 one obtains T oc (gN) 2 cx o~ 2 . 
We can now relate F to a rate of change of the width of 
the distribution a. In the case in which the expansion is 
assumed to be instantaneously diffusive, as done for ex- 
ample in to describe a lattice with random disorder, 
the diffusion rate is essentially T, and one obtains 



da 2 j dt oc a 



(8) 
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FIG. 6: Calculated mean coupling for doublets of states versus 
the distances between the states for a quasiperiodic lattice 
with A = 4J. The dominant couplings are those to states 
about one quasiperiod, d/(ki/ki — 1) ~ 4.4d, apart. 
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FIG. 7: Calculated probability of coupling for doublets of 
states in a quasiperiodic lattice with A = 3. 5 J (see text). 
The dashed and dotted lines are a guide to the eye for decays 
proportional to V-Eint and Ei n t, respectively. 



The large coupling of states separated by one quasiperiod 
shown in Fig. [5] suggests however that in our quasiperi- 
odic lattice the expansion might have a more coherent 
nature, i.e. the expansion is dominated by sequential 
hopping between states close in energy and separated by 
approximately one quasiperiod. In the case of such an 
expansion mechanism, the rate T represents instead the 
instantaneous velocity, and one obtains 

da/dt cx a~ 2 . (9) 

2) If instead Ei n t < SE, only a reduced number of mi- 
croscopic transitions are activated by the interaction, and 
the rate T has an additional dependence on Ei n t . For ex- 
ample, we numerically evaluated the behavior of the cou- 
plings between doublets of states for our quasiperiodic 
lattice with A = 3. 5 J (which corresponds to SE ss 3J), 
following the analysis presented in Refs. @, @|. To do 
this, we selected the strongest coupling \V- ikk \/AE ik for 
each initial state i and we calculated the probability to 
have such coupling active for a varying interaction en- 
ergy. We have then averaged this result over many initial 
states along the lattice. As shown in Fig. we find that 
the probability is about constant for Ei n t ^ SE, while 
it decreases at least as y/Ei n t for decreasing interaction 
energy. Since Ei nt oc a^ 1 , this implies that the related 
rate F oc (\V/ ikk | 2 ) can scale at least as a~ 3 . This changes 
the overall exponent in the equation of motion for a. For 
example, in the case of a coherent expansion one obtains 

da/dt cx a~ 3 . (10) 

This regime of slower subdiffusion is analogous to the one 
reported in Ref. for uncorrected disorder and dubbed 
"weak chaos" regime, in contrast to the opposite regime 
of " strong chaos" . 



FITTING THE AXIAL DYNAMICS IN THE 
EXPERIMENT AND IN THE NUMERICAL 
SIMULATIONS 

Let us now discuss in more detail the link between 
the microscopic model and the macroscopic expansion 
observed in the experiment and in the numerical simula- 
tions. We have derived above an equation of motion for 
the width of the form: da/dt — Ca~ p , where p is a posi- 
tive rational number and C is a constant which takes in 
account all the terms that do not depend directly on a. 
The integration of this equation with the initial condition 
tr(0) = do leads to 

<r(t)=<ro[l + -%ftt) , (11) 
V a(T o ) 

with a = This can also be expressed as 

a(t)=a (l + £\ , (12) 

where the parameter t represents an " activation time" . 
In general, three fitting parameters are therefore required 
to fit the data. For t >• t the initial size does not play 
a role in the expansion, two parameters in the fitting 
procedure are enough and the expansions shows a lin- 
ear behavior in the log-log scale. However, in order to 
decrease the uncertainty on the parameters and to ex- 
tract information also from the measurements at shorter 
times, we chose to fit also the data taken for t < to us- 
ing the three parameters fitting function. This gives also 
the possibility to confirm that the asymptotic regime has 
been reached by comparing to with the maximum obser- 
vation time. The characteristic value for the activation 
time to we fitted for the sets of data reported in the 
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paper is indeed around one second (except for the non- 
interacting case, where it is larger), i.e. one order of 
magnitude smaller than the maximum observation time. 
Note also that by fitting only the data points for t > t 
with the asymptotic form of Eq. [TT] we obtain exponents 
aii n which are in good agreement with the previous ones. 
One example of this comparison can be seen in Fig. |8ji. 
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FIG. 8: Typical time evolution of the axial width a for an in- 
teracting sample. The value for the diffusion exponent from 
a linear fit in log-log scale for asymptotic times au n (black 
dashed line) is in good agreement with the one from the three 
parameters fitting function a (red line). Panel a) showns 
the experimental data for Ei nt ~ 1.8; panel b) the numer- 
ical simulation for an interaction parameter /3 = 60 which 
corresponds to Emt ~ 1.8 

We also numerically solved the discrete ID DNLSE, 
Eq. 2 of Ref. 0] using the same parameters for the lat- 
tice as in the experiment. The interaction parameter /3 
of Ref. 0] is related to the interaction energy E int of this 
work by E int — {(3/2) J^i IV'il 2 ) where \tpi\ 2 is the proba- 
bility of finding an atom at the lattice site i. We report 
a typical curve obtained from a simulation up to 100 sec- 
onds in Fig. 4b. As already done with the experimental 
data, we extracted the exponent a by fitting the numer- 
ical data with the expression in Eq. QTJ We also checked 
that the exponent au n extracted by a fit of the data in 
the asymptotic regime alone is very close to a. 

Note that in the experiment we measure a ~ 0.06 even 
for Ei n t ~ 0. This is presumably due to a weak techni- 
cal noise on the quasiperiodic potential (e.g. vibrations 
in the retroreflecting mirror that creates the standing 
waves, or fluctuations of the laser frequencies). If we 
would assume that this value represents a constant bias 
for all measurements with interaction, then we would find 
experimental values for a that are closer to the theory. 



However, in literature there are no studies on the inter- 
play of noise and interaction in a disordered system that 
can fully justify this assumption. A study that we are 
currently performing on our system indicates that the 
diffusion exponents due noise and interactions roughly 
add up in the regimes explored in this work, but a con- 
clusive analysis is still missing. We therefore leave a more 
detailed discussion of this point to a future study. 



MODELING THE RADIAL DEGREES OF 
FREEDOM AND TEMPERATURE EVOLUTION 
IN THE EXPERIMENT 

The presence of the radial degrees of freedom in the 
experiment can now be taken into account. As shown 
in Fig. 1, the radial energy spacing Tiu r is about one or- 
der of magnitude smaller than the characteristic disorder 
energy SE. There are two main effects of the radial ex- 
citations to be considered: 

1) The 4-states overlap integral must be modified to 
include also the radial plane 



Iijkl = / f*f*fkfldx j iPriPrjiprkiPrld T 



(13) 



An increasing radial excitation of the system therefore 
results in a broadening of the radial wave functions, with 
a consequent decrease of the overlap integral and the as- 
sociated coupling term |Vyjy|. 

2) The possibility of populating excited radial modes 
softens the energy conservation requirement, since the 
radial level separation hoj r is smaller than the character- 
istic energy spacing of the quasiperiodic lattice SE (by 
approximately one order of magnitude in the present ex- 
periment). This implies that axial states that could not 
be coupled because of a too small E int in a ID system, 
might be at least partially coupled in presence of a radial 
excitation. 

We have not attempted so far a quantitative numerical 
study of this model including the radial modes. However, 
the considerations above can be used to derive qualita- 
tive expectations that can explain the experimental ob- 
servations. We can indeed conceive two possible limits 
for the radial temperature. The first one is the limit 
ksT r < huj r , which in our case (hui r /kB ~ 5 nK) would 
correspond to the minimum temperature achievable in 
the Bosc-Einstcin condensate. In this limit the excited 
radial population is negligible, and one has effectively a 
ID problem. However, this limit cannot be achieved as 
a stationary state, since the thermalization of the ax- 
ial and radial degrees of freedom will rapidly bring the 
radial temperature to a minimum value of the order of 
T r w SE/k b w 50 nK. The other limit is k B T r > SE, in 
which many radial states are excited. In this limit the 
energy mismatch that must be compensated by the cou- 
pling terms to provide the axial hopping is no longer of 
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order SE, but it becomes of order Tiu> r . The critical E int 
to reach the regime of strong chaos will therefore be hui r . 
The experiment is performed in an intermediate regime 
of temperature, ksT r SE, where we can expect that 
the critical energy for the strongly interacting regime lays 
somewhere in between SE and huj r . 
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FIG. 9: Typical time evolution of the radial temperature for 
a non-interacting sample (black squares) and an interacting 
one (blue triangles). 

This expectation is confirmed by the comparison be- 
tween the experimental and numerical energy depen- 
dence of the exponent a discussed in the paper. Note 
that the maximum temperature that can be reached in 
the experiment while preserving a single-band dynamics 
is of the order of 200nK. For higher temperatures we see 
a clear excitation of the second band of the quasiperiodic 
lattice, whose single-particle eigenstates are not localized 
for the range of parameters we explored. This excitation 
is detected as the appearance of rapidly moving tails in 
the density distribution. 

Finally, we discuss the evolution of the radial temper- 
ature T r in the experiment for samples that are initially 



prepared as Bose-Einstein condensates, i.e. at the low- 
est measurable temperatures, T r < 50nK. The typical 
evolution is reported in Fig. |H1 In absence of interac- 
tion we see a slow heating with a rate of the order of 1-3 
nK/s, presumably due to pointing or amplitude noise on 
the lattice beams. The interacting samples show instead 
a larger initial T r s» SE/ks, which presumably arises 
from a thermalization of the axial and radial degrees of 
freedom. The following heating during the full evolu- 
tion time is typically 3-5 nK/s, hence larger than the one 
of non-interacting samples. This excess heating presum- 
ably arises from the same axial noise on the quasiperiodic 
lattice that causes the slow expansion of non-interacting 
samples, see Fig. 1 in the paper. This noise might be due 
to variations in the laser wavelengths and to vibrations 
of the retro-reflecting mirror that creates the standing 
waves, which result in a spatial movement of the two lat- 
tices and a consequent excitation of the axial dynamics 
even in the case of orthogonal localized states. In pres- 
ence of interaction the axial excitation is transferred to 
the radial degrees of freedom, resulting in an increased 
heating. 
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